Self-assembly of amphiphilic Janus particles at planar walls: A density functional 

study 
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We investigate tiie structure formation of amphiphilic molecules at planar walls using density 
functional theory. The molecules are modeled as (hard) spheres composed of a hydrophilic and hy- 
drophobic part. The orientation of the resulting Janus-particles is described as a vector representing 
an internal degree of freedom. Our density functional approach involves Fundamental Measure The- 
ory combined with a mean-field approximation for the anisotropic interaction. Considering neutral, 
hydrophilic and hydrophobic walls, we study the adsorption of the particles, focussing on the com- 
petition between the surface field and interaction-induced ordering phenomena. Finally, we consider 
systems confined between two planar walls. It is shown that the anisotropic Janus interaction yields 
pronounced frustration effects at low temperatures. 



I. INTRODUCTION 

Amphiphilic molecules are generally composed of a po- 
lar, hydrophilic (i.e. " water- loving" ) head and at least 
one hydrophobic (i.e., "fat-loving") tail, typically a car- 
bon chain. The presence of these two ingredients, in 
combination with an aqueous solvent, yields a large vari- 
ety of self-assembled structures on different length scales, 
ranging from molecular-size micelles to mesoscopic mem- 
branes, bicontinuous foams and lamellar phases. From an 
application point of view, amphiphiles are used in a va- 
riety of contexts, e.g., to reduce the surface tension in 
complex mixtures (such as water and oil), as templates 
to create nanoporous materials^ and, more recently, to 
stabilize bundle and network formation in solutions of 
carbon nanotubes^'^. 

From the theoretical side, the self-assembly of am- 
phiphilic molecules has been investigated by a variety of 
approaches and modelsl'^~, including lattice gas systems^, 
Ginzburg-Landau theory'^, density functional studies of 
entropic models (where the surfactant is represented by 
a sphere plus an infinitely thin rod^^— , or by a dimer-ii) , 
and off-lattice simulation o^^d^ of flexible bead-spring 
moleculesii, with and without explicit solvent. A par- 
ticularly simple, implicit-solvent model has been pro- 
posed by Tarazona et al. already in 1995^^. Within this 
model, the amphiphilic particle is a sphere composed of 
two hemispheres, one being hydrophilic and the other 
one hydrophobic. The solvent-mediated interaction is 
then taken into account via an effective directional po- 
tential involving the orientation vectors defined by the 
symmetry axis of the spheres. Despite these simplifica- 
tions (which completely neglect geometrical factors such 
as relative size of the head group etc.), the model is ca- 
pable of describing bilayer-, vesicle-, and micelle forma- 
tion. From todays perspective, one would call this sim- 
plified amphiphilic molecules rather as "Janus spheres". 
This term was originally proposed by Casagrande et al. 
in 198&i^ to describe spherical glass particles with one 
of the hemispheres being hydrophilic and the other one 
being hydrophobic. Since then, and as anticipated by 



P. G. de Gennes on the occasion, of his Nobel lectureii, 
the surface properties of these Janus particles have be- 
come an area of great interest on their own. Gener- 
ally, Janus particles are composed of at least two phys- 
ically (or chemically) distinctive surfaces, and there has 
been significant progress to synthesize such particleai^"— 
(whose interactions can also be, e.g., dipolar or magnetic 
in character—"—). Theoretically there are several models 
"on the market" which are similar in character to that 
proposed by Tarazona et al}^ , but have been indepen- 
dently proposed specifically for Janus particle a^^i^^'^'' . 
Furthermore, in a recent study of Sciortino et al^ Janus 
particles are described within a "patchy-particle" model. 

In the present paper we are interested in the self- 
assembly of amphiphilic Janus particles at planar in- 
terfaces and in strongly confined situations (slit-pores). 
Indeed, within the more general context of surfactant 
assembly, the impact of surfaces is important in a va- 
riety of applications, including synthesis of nanoporous 
materials, thin film deposition for lithographic processes 
and devices^, and enhancement of chemical reactionSi^. 
An additional attractive factor of studying systems with 
surfaces is that the self-assembled structured can be 
experimentally studied using e.g., neutron reflectom- 
etry (targeting the thickness of the layer)—, grazing 
incidence small-angle neutron scattering'^^ and atomic 
force microscopy (targeting the lateral structure) and 
by studying adsorption isotherms--. Moreover, there 
is a strong fundamental interest to understand the 
self-assembly process of Janus-like particles at inter- 
faces as a bottom-up process for the design of future 
nanomateriala^. 

Contrary to the extensively discussed case of true sur- 
factant molecules at surfaces, the self-assembly of Janus 
particles at surfaces has so far rarely been considered, 
an exception being a recent study of Hirose et al.^^ who 
used a macroscopic theory based on the Young's equa- 
tion. In the present paper we use a microscopic approach, 
that is, classical density functional theory-''!'^*. Follow- 
ing the original work of Tarazonapi^, we employ a mean- 
field approximation to treat the anisotropic part of the 
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interactions, whereas the repulsive (hard sphere) contri- 
bution is treated on a more sophisticated level. Specifi- 
cally, we employ the so-called Fundamental Measure the- 
ory (FMT)22 which has turned out to be extremely suc- 
cessful for the description of inhomogeneous hard sphere 
systems^i^. As a starting point, our investigations focus 
on planar structures such as bilayers. The key question 
is to determine thermodynamics conditions under which 
self-assembly arises at the surface, as compared to the 
corresponding bulk system. Moreover, we explore the 
impact of different surface properties concerning, in par- 
ticular, their hydrophilic or hydrophobic character. 

The remainder of the paper is organized as follows. 
In Sec. Ill Al we introduce the model and also discuss its 
main features as compared to other recently proposed 
models of Janus-like systems. The density functional for- 
malism including the FMT contribution is presented in 
Sec. IIIBi accompanied by an appendix which describes 
relevant technical issues. Numerical results are presented 
in Sec. lIIIi where we focus on the case of single walls, but 
briefly discuss also the case of Janus particles in slit-pore 
confinement. We close the paper with a summary in 
SecIVl 
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FIG. 1: Left: Illustration of typical pair configurations within 
our model. Right; the first column of the table gives the ori- 
entational contributions to the pair energies [i.e., the term 
(ui — U2) ■ f2i] according to Eq. 1)2. 3p . The second and third 
columns refer to the possible extensions of the model accord- 
ing to Eqs. (|2.6p and (|2.7p . respectively. 



II. THEORY 



where 



A. Model 

In this study we employ a simple, coarse-grained model 
of an amphiphilic system originally suggested by Tara- 
zona and coworkersi^. In this model the amphiphilic 
molecules are represented by spherical particles consist- 
ing of two hemispheres, one being hydrophilic and the 
other one being hydrophobic. Thus, the "molecules" 
rather resemble Janus particles. The solvent, which is 
omnipresent in real amphiphilic solutions, is treated im- 
plicitly. The total pair potential between two Janus-like 
particles at positions ri and Vi subdivides into a hard- 
sphere (HS) and an anisotropic contribution. 



i(ri2,ui,ii2) 



-HS 



(»-12) + (/-^ (ri2,Ui,U2) , (2.1) 



where ri2 = ri — r2, and the HS potential is given by 



</>«^ (ri2) 



00, ri2 < cr, 
0, ri2 > cr. 



(2.2) 



with ri2 = |ri2| and cr being the HS diameter. Further, 
Ui and U2 arc unit vectors denoting the orientation of 
the sphere. Specifically, (i — 1,2) points from the 
hydrophobic to the hydrophilic side of particle i. The 
vectors u = u (w) are parameterized by u; = ((/j, 0), with 
(/? = [0,27r] and 9 = [0,7r]. The anisotropic interaction is 
defined as 



6^ (ri2, ill, U2) = 01 (ri2) (£ii - U2) • r2i. 



(2.3) 



|^6exp(-A(ri2 -cr))/ri2, ri2 > cr. 

In Eq. (|2.4p the parameters C and A measure the cou- 
pling strength and the (inverse) range of the interaction, 
respectively. If not stated otherwise, we set Acr = 3. 

A characteristic feature of the anisotropic potential 0^ 
is that it lacks of any coupling between the orientations \\i 
of the two involved particles (i = 1, 2). In fact, Eq. ()2.3p 
represents only the lowest-order anisotropic contribution 
of a more general expansion of the pair potential be- 
tween linear molecules into spherical harmonics-S (see 
also Eq. (|2?5|) below). 

In Fig. [1] we sketch the most important configurations 
of two particles; the corresponding pair energies accord- 
ing to the angle-dependent terms in Eq. (12. 3p are given 
in the related table (first column). Within our model 
the particles prefer to be orientated in opposite direc- 
tions, such that the hydrophobic sides point towards one 
another. The opposite configuration (with facing hy- 
drophilic sides) is the energetically most unfavorable one. 
Further, parallel orientations are energetically neutral. 
The strong energetic preference (or penalty) of configura- 
tions with facing hydrophobic (hydrophilic) sides mimics 
the effects expected in a real system which would include 
water as a solvent. Clearly, the water molecules will pref- 
erentially adsorb at the hydrophilic side of each particle. 
The resulting steric exclusion yields an effective repul- 
sion of the hydrophilic sides of neighboring Janus parti- 
cles. On the other hand, the fact that the hydrophobic 
sides dislike water effectively favors configurations where 
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these sides point towards each other, such that the con- 
tact with water is minimized. An additional source of 
attraction between hydrophobic sides can arise from a 
functionalization with, e.g., Thiol. An example has been 
described in Ref. HI where Thiol (and other molecular 
groups) is used to functionalize the gold-coated side of a 
polystyrene sphere. The functionalization yields a deple- 
tion area with respect to water and thus, an additional 
effective attraction between the hydrophobic (gold) hemi- 
spheres. 

It is interesting to briefly consider the subsequent 
terms [i.e., beyond Eq. (j2.3l) ] in the general expansion 
of the anisotropic pair potential. These terms readi^ 

(j)^ (ri2, ui, U2) = 01 (ri2) (ui - U2) • f2i 

+9102 (7-12) Ql (f21, Ui, U2)+q203 {ri2) Q2 (f21, Ui, U2) , 

(2.5) 

where 

Ql (f21, Ui, U2) = ( - (Ul ■ f2l) (U2 • f2l) 



1 - (Ul • f2l)^ V 1 - (U2 • f2l)^C0S {lPi - ip2) 



The remaining parts of the spheres only induce steric re- 
pulsion. This picture is quite different from our model, 
which should rather be compared to a particle with two 
patches, one mimicking the hydrophilic and the other 
one the hydrophobic part. Moreover, in our model the 
strongest attractive (repulsive) interactions occur when 
the particles are aligned in an antiparallel (parallel) man- 
ner. The one-patch model in Ref. 2& concentrates on the 
attraction in "antiparallel" configurations, i.e. configu- 
rations where the two patches become coupled through 
the square-well zone. 

Our main goal in the present paper is to study the 
amphiphilic system in the presence of a planar surface 
(located at z = 0). The simplest surface model is a hard 
wall, yielding the surface potential 



0™''f (z) = 



00, z < (t/2, 
0, Z> cr/2. 



(2.8) 



To include effects of preferential adsorption we also con- 
sider potentials of the form 



-Ul • U2, (2.6) i,surf 



>,u) 



and 



(32 (f21,Ui,U2) - (3(U1 • r2l) +3(u2-f2l) " 

(2.7) 

In Eq. (|2.5|) . 02 (''12) and 03 (ri2) are distance-dependent 
functions that could be set, for example, equal to the 
Yukawa-like function 0i (ri2) defined in Eq. (12. 4p . The 
implications of the additional, angle-dependent contri- 
butions are illustrated in the center and the right col- 
umn of the table in Fig. [T] Equation (j2.6p adds a term 
(Ql) involving the negative scalar product of the orien- 
tation vectors; this term clearly prefers parallel orien- 
tations irrespective of the configuration of the particles 
relative to the connecting vector. The term Q2 defined in 
Eq. (|2.7p favors parallel orientations perpendicular to the 
connecting axis. Here we focus on the lowest-order term 
[cf. Eq. (|2.3I) ] which describes the main characteristics of 
our amphiphilic system^-. 

To complete the discussion on the interaction between 
the Janus particles we briefly compare the present model 
with other models that were recently proposed in the 
literature. The potential suggested by Hess et al^ in- 
volves all types of angle-dependent terms appearing in 
Eqs. (|2.3p - (|2'7p . However, the parameters in Ref. [2^ were 
chosen such that parallel orientation of the Janus spheres 
is preferred, contrary to our model favoring antiparallel 
alignment (see Fig. [ij . Another Janus potential has been 
suggested by Sciortino et al?^ . Within this model, the 
(spherical) particles have one patch; the patches of neigh- 
boring particles then interact via a square- well potential. 



00, z < <j/2, 

C2 exp (-A2 (2 - cr/2)) u • 62/2, 2: > 0-/2, 

(2.9) 

where the unit vector points along the z-direction. 
The sign of the parameter C2 determines the preferred 
orientations of the particles. Specifically, for negative 
(positive) values of C2 the surface prefers the hydropho- 
bic (hydrophilic) side of the Janus particle. Finally, the 
range of the surface potential is controlled by the param- 
eter A2. For simplicity, we set A2 = A = 3a~^. 



B. Density functional theory 



The central quantity in our study is the normalized 
singlet density 



p (r, u:)^p (r) a (r, a;) , / (fuja (r, a;) = 1, (2.10) 



which depends on both, the position and the orientation 
of the particles. In the equation above, cPlo — dipdOsind. 
Following other DFT studies of molecular fluids (see, e.g., 
Refs.— i^i^) we write the singlet density as a product of 
a position-dependent number density p (r) and an orien- 
tational distribution function a(r,u}). The equilibrium 
singlet density is obtained from the grand-canonical den- 
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sity functional f2 [p] , which is given as 



Q [p, a] = /3-1 J d?Tp (r) (in {p (r) A^) - l) 

+ /3"^ j j (fr(fujp (r) a (r, lj) In (47ra (r, lj)) 
+ jj dPrdPijjp{v)a{v,u})(lf'''^ (r,u(a;)) 



+ d^rd^r'Sujp (r) p (r') 
a(r,a;)^i(r-r')^^— ^-uCo;). (2.11) 



The first two contributions on the right side of Eq. (|2.1ip 
represent the ideal gas part of the free energy, with 
P = l/k^T (where fee is the Boltzmann constant, T the 
temperature) and A is the thermal de Broglie wavelength. 
Further, the third and the fourth term include contribu- 
tions from the external (surface) potential and the chem- 
ical potential, respectively. The fifth term arises from the 
HS interactions, which we treat by FMT— as described 
in more detail below. Finally, the last term in Eq. (j2.1ip 
involves the anisotropic interactions (|2.3p . We treat this 
free-energy contribution by a mean-field approximation 
(which can be rationalized from a A-expansion^ of the 
pair potential). In this framework, the pair correlation 
function is set to one for separations ri2 > a. Thus, cor- 
relation effects beyond the core condition are neglected 
within the anisotropic contribution to the density func- 
tional. 

Regarding the HS free energy, we use the White-Bear 
mark II version^ of FMT. In this framework the function 
^WBii appearing in the fifth term on the right side of 
Eq. (|2.11l) is given by 



(/)«^^^({n„ (r)}) = -no (r)ln(l-n3 (r)) 



1 



V'3 ("3 (r)) 



^ ni (r) n2 (r) - ni (r) ■ n2 (r) 

1 - 713 (r) 
^ 712^ (r) - 3712 (r) na^ (r) 
24^(1-713 (r))' 

In Eq. (EH]), 

V'2 {ns (r)) = ^^(2*^3 (r) - (r) 



'<p2 (?^3 (r)) 



(2.12) 



and 

ip3 ("3 (r)) 



m (r) 

+ 2(l-773(r))ln(l-773(r))), (2.13) 



(2773 (r) - 3773" (r) + 2n3^ (r) 



7732 (r) 



2(l- 773 (r))'ln(l- 773 (r))). (2.14) 



In the expressions above, the scalar functions ria (r) with 
a = 0, 1, 2, 3 and the vector functions (r) with 7 = 1,2 
are weighted densities defined as 



77a (r) 
Hq (r) 



= / d\'p (r - r') 



UJa (r') 
u^a (rO 



(2.15) 



The weight functions are composed of Dirac delta dis- 
tributions 6{x) and Heaviside step functions Q{x). Ex- 
plicitely, one has 



c:;3 (r) = e(i?-|r|),cS2 (r) = 6 {R - \r\) 



ui (r) 



UJ2 (r) 
AttR 



, ^0 (r) 



^2 (r) 
47ri?2 ' 



(2.16) 



(2.17) 



(2.18) 



where R = (j/2. 

The equilibrium density is found by minimization^^ of 
the grand canonical functional given in Eq. (j2.1ip . that 
is 



dn [p, < 



Sp{r) 



Po(r) 



sn [p, 



Sa (r, iv) 



0, 



= 0. 



(2.19) 



(2.20) 



Due to the mean-field approximation of the anisotropic 
interaction in the density functional [see last term in 
Eq. ()2.1ip ]. combined with the fact that this interaction 
is linear in the orientations [see Eq. ()2.3p ]. minimization 
with respect to a (r, a;) yields the simple, explicit expres- 



ao (r,a;) 



exp (a (r) • u (a;)) 
/ d^w'exp (a (r) • u (a;')) 



(2.21) 



where 



a (n) = /3 y dV2p (r2) 01 (ri2) f 12. (2.22) 

Contrary to the orientational distribution, the 
position-dependent number density has to be calculated 
numerically. Technical details are given in the Appendix, 
where we also give the exact boundary conditions by 
which the numerical results can be checked. To restrict 
the numerical effort we focus on density distributions 
that depend only on z. Clearly, this strategy implies 
the possibility of missing other, and maybe energetically 
preferable solutions. Indeed, by considering both, z- 
dependent and r-dependent density distributions, Tara- 
zona et al^ showed that the present model yields planar 
structures ("membranes"), but also spherical structures 
("vesicle" and "micelles") under bulk conditions. How- 
ever, in the presence of surfaces, a full three-dimensional 
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calculation would demand considerable computational 
cost. Therefore, we focus on density distributions with 
the symmetry suggested by the planar wall and assume 
translational symmetry in the other spatial dimensions. 
As a consequence, the orientational distribution also sim- 
plifies to a (r, <jj) — a{z, 9). 

Having this in mind, it is useful to introduce various 
^-dependent orientational order parameters. Specifically, 
we define 

/i(z) = 27r/" d6lsin6'Q!(z,6')cos6' (2.23) 
Jo 

and 

/i2 (z) = TT / d9shi9a{z, 9) (Scos^ei - 1) . (2.24) 
Jo 

The function h{z) describes the "polarization". It can 
take values between — 1 and 1 corresponding to complete 
alignment in negative (positive) z-direction at the posi- 
tion z considered. The latter, second-rank order param- 
eter h2{z) also describes the ordering along the z-axis, 
but it is independent of the direction. Within our (mean- 
field) approximation, the function h2{z) is strongly cou- 
pled to h{z), that is, h{z) ^ automatically yields 
/i2 {z) ^ and vice versa (see also Fig. H]). Hence, we 
mainly consider the z-dependent "polarization", h{z). 



III. RESULTS AND DISCUSSION 

A. Bulk behavior 

As a starting point we briefly consider the bulk sys- 
tem of Janus spheres, which may be characterized by 
the packing fraction, -q = {Tr/6)pa^, and the reduced 
temperature T* = k-eT/C [with C being the coupling 
constant appearing in Eq. (|2.4p ]. With the latter def- 
inition, the limit T* — )■ oo corresponds to a pure HS 
system. Therefore, at high T* and densities 77 < 0.45, 
minimization of the density functional (j2.11l) just yields a 
homogeneous, isotropic phase characterized by constant 
number density, p{z) — p, and constant orientational dis- 
tribution, a{z^9) ~ l/47r. Upon lowering the tempera- 
ture, inhomogeneous, "membrane" -like solutions can ap- 
pear. By "membrane" we refer to an isolated bilayer of 
Janus particles, which is infinitely extended along the x- 
and j/-directions, and where the two sub-layers are ori- 
ented towards each other in the energetically preferable 
configuration (with facing hydrophobic sides). The ap- 
pearance of such membranes was already shown in the 
original paper of Tarazona et ali^. To detect the cor- 
responding temperature range in the present study, we 
start the numerical minimization with an initial guess for 
p{z) composed of two Gaussians. At high temperatures, 
this initial profile quickly disappears during the itera- 
tion procedure. For low values of T*, on the other hand. 






0.1 0.2 0.3 

FIG. 2: Bulk phase diagram. The solid line indicates the 
estimated transition temperatures (see main text) between 
the homogeneous, isotropic high-temperature phase and an 
ordered low-temperature phase characterized by membranes. 
The dashed line corresponds to the vapor-liquid coexistence 
curve resulting from a density function calculation with a 
MMFT treatment of the anisotropic interactions (yet under 
the assumption of a homogeneous isotropic state) . 



the procedure yields indeed formation of periodically re- 
peated bilayer structures, characterized by a sharp de- 
cay of the density outside of the bilayers, and a depletion 
zone between the membranes. The temperatures separat- 
ing both regions are indicated by the solid line in Fig. [2] 
However, the precise location of this line has to be con- 
sidered with some caution since the inhomogeneous low- 
temperature profile did not converge to a truly stable 
result. Therefore, we should consider the line in Fig. [2] 
rather as an upper limit of the temperature range related 
to membrane formation. We also note that our results 
differ from those obtained in Ref.fisl where the transition 
temperatures were found to be higher (the focus ir^i^ was 
on lower densities, though). These deviations might be 
due to the fact that the HS part of the density functional 
used in Ref. [l^ was treated in weighted-density approxi- 
mation, rather than by FMT as in the present study. 

We also note that the present functional does not 
predict a gas-liquid transition of the Janus-particle sys- 
tem, as it was found, e.g., in a recent simulation study 
of a related model system consisting of patchy Janus 
particles^-. Within our study, the absence of condensa- 
tion is a consequence of the mean-field treatment of the 
anisotropic interactions; indeed, the corresponding term 
in the density functional [see last term on the right side of 
Eq. (|2.11l) ] is zero in a fully disordered system. Interest- 
ingly, if we focus on such a homogeneous, isotropic system 
and use a modified mean field approximation (MMFT) 
for the anisotropic part (for applications of the MMFT 
to related systems, see ReisM^^), we do find a gas-liquid 
phase transition, indicating that the Janus interaction is 
effectively attractive. The corresponding phase coexis- 
tence points are indicated by the dashed line in Fig. [21 
Here we do not explore this aspect further, the main rea- 
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FIG. 3: (a) Density profiles and (b) order parameters at 77 = 
0.05 and various temperatures. Here and in the subsequent 
figures, the hard wall is located at 2; = 0, and the range 
parameter is set to Act = 3. 



son being that the use of MMFT in the presence of sur- 
faces (which are the focus of our study) would have im- 
plied a drastic increase of numerical effort. In the follow- 
ing discussion of surface effects (on the basis of the pure 
mean-field approximation) we therefore consider temper- 
atures and packing fractions in the whole range above the 
solid line in Fig. [21 i.e., where the bulk system is disor- 
dered. 



B. Surface effects at neutral walls 
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FIG. 4: Density profiles (a) and order parameter functions 
h{z) (b) and h2{z) (c) at r; = 0.3 and various temperatures. 



We start the discussion of surface effects with the case 
of hard ("neutral") walls [see Eq. (H^H])]. In Fig. ISJa) we 
present density profiles for a dilute system characterized 
by r/ = 0.05 and different temperatures T* > T^, where 
~ 0.079 is the temperature related to membrane for- 
mation in the bulk (see Fig. [5]) . The corresponding or- 
der parameter function h{z) [sec Eq. (I2.23p ] is plotted in 
Fig. [3Jb). At T* — 10, the density profile agrees within 
the numerical accuracy with that of a HS fiuid, consistent 
with the negligible values of the order parameter. By de- 
creasing the temperature, configurations with neighbor- 
ing hydrophobic sides of the particles become more and 
more favorable, yielding orientational ordering of the par- 
ticles as visible in Fig.[31[b). At the same time, the density 
profiles change from the typical behavior of a dilute HS 
system, where the density maximum occurs directly at 
the wall, towards a soft, loosely packed structure where 
most of the particles agglomerate somewhat away from 
the waff [see, e.g., data at T* = 0.08 in Fig.lHfa)]. More- 
over, the distance between the first two layers indicated 
by the two close maxima in p(z) is smaller than one parti- 
cle diameter. Analyzing the corresponding orientational 
order parameters we find that, in the first layer (located 
at z w la) the particles are oriented such that their hy- 
drophilic side points preferentially to the wall, whereas 
in the second layer (at z 1.8 — 2cr), they tend to ori- 
ent in the opposite way. Thus, despite the low density 



considered, a significant degree of orientational ordering 
is already present. 

As expected, both the translational and the orienta- 
tional ordering becomes more pronounced upon an in- 
crease of the density. This is illustrated in Fig. |4] where 
we plot density and orientational profiles for the exem- 
plary case 77 — 0.3. At the lowest temperature considered, 
T* = 0.28, the Janus particles close to the wall are almost 
completely arranged into a double layer, as reflected by 
the high values of the first, pronounced maxima in p{z) 
located at 2: = 0.5a and z = 1.5a. The double-layer 
formation is accompanied by a high degree of "polar" 
orientational order, as seen from the pronounced change 
of the function h{z) [see Fig. IH^b)] from negative values 
at z w 0.5(7 (indicating that the hydrophilic side tend to 
point towards the wall) to positive values at z w 1.5a. 
Part (c) of Fig. m additionally shows the function h2{z) 
defined in Eq. (|2.24p : however, as anticipated in Sec. Ill Bl 
this function just follows the oscillations of h{z) and will 
thus not be considered further. Combing back to the den- 
sity profile p{z) [see Fig. IH^a)] we see that, at distances 
z beyond the bilayer at the wall, the third maximum ap- 
pears only at 2 w 2.7a, indicating that the next layer 
is slightly shifted towards larger separations. This is a 
result of the repulsion between the hydrophilic sides of a 
particle in the second, and one in the third layer. Indeed, 
the preferred orientation of the third layer becomes ap- 
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parent from the positive value of h{z) in Fig. SJ^b). We 
stress that aU of these orientational ordering effects occur 
at temperatures where the corresponding bulk system is 
still homogeneous (see Fig. [5]). 

Within the framework of neutral walls, we have also 
considered the influence of the range parameter A within 
the Janus-particle interaction [see Eq. (|2.4p ]. It turns 
out, however, that the latter has no crucial effect for 
the systems considered here. Upon increase of A (i.e., 
decrease of the interaction range) one merely observes a 
decrease of the density maxima close to the walls, and a 
damping of the oscillations at larger distances. 



C. Influence of the surface fields 

We next consider the ordering behavior in systems, 
where the pure confinement effect induced by a neutral 
wall is supplemented by surface fields preferring (locally) 
a specific orientation of the Janus particles [see Eq. (j2.9|) ] . 
We start with a wall preferring the hydrophilic sides (such 
as silicaSi). This situation implies that the surface field 
supports the orientation already found close to neutral 
walls (see Sec. IIII Bl) . Exemplary density profiles and 
order parameters are shown in Fig. [51 where the bulk 
density i] = 0.3, and the fluid-fiuid and fluid-surface in- 
teractions (as well as the corresponding range parame- 
ters) are of the same magnitude. Due to the hydrophilic 
(and thus, supportive) character of the wall, the density 
in the contact layer, as well as the corresponding orien- 
tational order parameter, is even enhanced as compared 
to the corresponding system with neutral walls. This can 
be seen when we compare, e.g., results for the tempera- 
ture T* = 0.31 in Figs.[5][a) and (b) with corresponding 
ones at T* — 0.3 in Fig. 2] Contrary the neutral- wall 
case, however, the second density maximum at the hy- 
drophilic wall is much smaller than the first. We under- 
stand this behavior as a consequence of the fact that, in 
the second layer, the hydrophilic orientation preferred by 
the wall competes with that dictated by the fluid-fluid- 
interaction; the latter rather favors the hydrophobic side 
orienting towards the wall. The competition is also re- 
flected by the asymmetric shape of the maximum in the 
function h{z). Finally, considering the third density max- 
imum and comparing with Fig. |4l we flnd that this (and 
the subsequent) peak(s) is higher and that the depletion 
effect (relative to the second layer) is less pronounced 
than at a neutral wall. On the other hand, the density in 
the depletion zone is lower in the case of the hydrophilic 
(than at a neutral) wall. We conclude that, despite its 
short-range character [see Eq. (|2.9p ]. the surface field is 
still effective even at fairly large separations from the 
wall. 

Completely different behavior is observed at a surface 
preferring the hydrophobic side of the particles. This sit- 
uation is depicted in Fig. |5] where 77 = 0.3 and we have 
chosen, for the purpose of illustration, a fluid-wall cou- 
pling parameter twice as large as that of the fluid-fluid 




FIG. 5: (a) Density profiles and (b) order parameters at ?7 = 
0.3 and a hydrophilic wall (C2 = C, X2cr = Act = 3). The 
inset in part a) shows the density profile on an enlarged scale. 



interaction [i.e., C2 = — 2C in Eq. (|2.9p ]. As a result 
of the dominant surface field, the bilayer formed at neu- 
tral walls (see Fig. S]) completely breaks down, and one 
rather observes formation of a monolayer at sufficiently 
low temperatures. Moreover, this monolayer is build by 
particles whose hydrophilic side points away from the 
wall, contrary to what has been observed before. The 
subsequent layer then has the reverse orientation and is 
shifted to slightly larger distances (z w 1.65) as a result 
of the Janus repulsion between first and second layer. 
Behind these first two layers, one observes the typical os- 
cillatory density and orientation profiles, indicating the 
tendency for bilayer formation consisting of oppositely 
oriented Janus particles. 

Clearly, all of these effects depend on the ratio between 
the fiuid-wall and fluid- fluid coupling parameter, C2/C, 
and on the range parameter A2. Specifically, upon in- 
crease of IC2/CI one observes an increase of the extrema 
of both, p{z) and h{z), and an enhancement of the de- 
pletion areas. Similar effects emerge upon an increase of 
A2, that is, an decrease of the range of the surface field. 
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Finally, we consider systems confined between two pla- 
nar walls, that is, in a slit-pore geometry. For simplicity, 
we focus on the case of neutral walls. For large wall sepa- 
rations Lz, one expects the structure at either surface to 
become decoupled from that at the other surface, yield- 
ing bulk like-behavior (i.e., p(z) = Pbuik, h{z) = 0) in be- 
tween the two walls. Indeed, we have explicitely checked 
that under such weakly confined conditions, the single- 
wall behavior discussed in Sec. IIII Bj is recovered at each 
of the walls. For smaller L^, pronounced confinement ef- 
fects appear. Exemplary density and orientation profiles 
are plotted in Fig. [71 where = 7cr (and 77 = 0.3). At 
high temperatures (i.e., close to the HS limit) the system 
is almost bulk-like in the center of the slit-pore. Decreas- 
ing the temperature we observe layer formation (which 
is typical for any confined system), accompanied by the 
development of orientational order particularly close to 
the walls. At the wall separation considered, the "polar- 
ization" at low temperatures is asymmetric in the sense 
that particles in the left contact layer point antiparal- 
lel to those in the right contact layer. This is consistent 
with the fact that the low-temperature system consists of 
seven layers, a structure which allows for three full bilayer 



FIG. 8: (a) Density profiles and (b) order parameters at 77 = 
0.3 and various temperatures. The walls located at z = Oct 
and z = 3.5a. 

structures (composed of oppositely oriented particles as 
discussed in Sec. IIII B[) plus one single layer. 

Two more examples are shown in Figs. [S] and [HI corre- 
sponding to the cases = 3.5a and = Sti, respec- 
tively. In the first case, the high-temperature system is 
characterized by three layers of particles, with the mid- 
dle layer being rather thick [see Fig.[5l^a)]. Upon lowering 
the temperature, the anisotropic fluid-fluid interactions 
yields a splitting of the middle peak, reflecting that par- 
ticles in the middle layer tend to arrange in a "buck- 
led" structure where neighboring particles are somewhat 
shifted to each other with respect to the z-direction. The 
corresponding orientation profile in Fig.|Sl[b) reveals that, 
in this buckled middle layer, the particles arrange in an 
antiparallel way. Note that the resulting in-plane ar- 
rangement is not particularly unfavorable, since accord- 
ing to our model, side-by-side configurations are energet- 
ically "neutral" (see Fig. [T]). By assuming this rather 
complex structure the system at = 3.5a overcomes 
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FIG. 9: Same as Fig. [S] but for walls located at z = Oct and 
z = 3ct. 
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FIG. 10: Normal pressure as function of the wall separation 
for two temperatures. 



frustration effects. Even stronger confinement, as it is the 
case at — 3cr (see Fig.lHl), then yields three pronounced 
layers of particles, with the contact layers pointing in op- 
posite direction. However, contrary to the situation at 
Lz = 3.5(7, the order parameter directly at the position 
of the middle layer is zero. Only slightly left or right 
of the center one finds a preferred Janus-like orientation. 
We therefore regard this case as a frustrated system. 

The different microscopic configurations appearing in 
the strongly confined systems in dependence of Lz give 
rise to pronounced oscillations of the normal pressure Pz 
[as calculated from the contact theorem, see Eq. (lASp j. 
The importance of this quantity stems from the fact that 
it is experimentally accessible, e.g., by colloidal-probe 
atomic force microscopy^. Results for P^{L^ at two 
temperatures are shown in Fig. 1101 As it is typical for 
confined, dense systems of spherical particles, the oscilla- 
tions have a period of about one particle diameter. Com- 
pared to the high-temperature situation, the anisotropic 
interactions not only strongly enhance the amplitude of 
the oscillations; they also lead to slightly asymmetric 
peak shapes and to a shift of the oscillations. Indeed, at 



low temperatures, the maxima of Pz(iz) occur at mul- 
tiples of the particle diameter, consistent with the pre- 
viously discussed frustration effects, e.g., at Lz = 3cr. 
On the other hand, the complex structure seen, e.g., at 
iz — 3.5(7, corresponds to a minimum of the normal 
pressure curve. 



IV. SUMMARY 

In this paper we have used classical density functional 
theory to explore the structure formation of amphiphilic 
Janus particles at planar surfaces. Our density functional 
approach involves a sophisticated (FMT) treatment of 
the repulsive (HS) interactions, whereas the anisotropic 
interactions are treated on a mean-field level. One key 
finding of our study is that, due to the presence of a 
surface, significant translational and orientational order- 
ing related to bilayer formation occurs under conditions 
where the bulk system is still homogeneous and isotropic. 
Thus, the surfaces seem to strongly support the structure 
formation, even when this surface is just a neutral (hard) 
wall. Moreover, we have shown that the details of the in- 
homogeneous structure at the wall can be "tuned" by 
varying the surface potential. Indeed, walls preferring 
the hydrophilic part tend to enhance the bilayer struc- 
ture seen already at neutral walls, whereas hydrophobic 
walls typically induce a competition between mono- and 
bilayer structures. We note that the degree of hydropho- 
bicity can be experimentally tuned, e.g., by coating sil- 
icon wafers with polymer films of varying thickness'*-. 
Finally, we have considered confinement effects emerg- 
ing from the presence of two planar (neutral) surfaces. It 
turns out that, for specific wall separations, there are pro- 
nounced frustration effects stemming from the interplay 
between the fiuid-surface potential, which prefers planar 
layering, and the fiuid-fluid potential preferring bilayers 
with depletion areas in between. This competition is also 
visible in an experimentally accessible quantity, that is, 
the normal pressure as function of wall separation. 

Clearly, the present study is only a starting point for 
a more systematic investigation of the impact of surfaces 
on the self-assembly of amphiphilic Janus particles. From 
a physical point of view, one main drawback of our cal- 
culations is the restriction to planar, self-assembled ge- 
ometries. Especially at low densities one would also ex- 
pect the occurrence of spherical structures such as vesi- 
cles and micelles; the most "stable" structure could then 
be selected by comparison of the related free energies. 
Furthermore, in experiments of elongated amphiphilic 
molecules (rather than amphiphilic spheres) at surfaces, 
both planar and spherical structures are observed^^^, 
suggesting that different structures could also occur for 
the spherical (Janus) case. The strategy to include spher- 
ical self-assembled structures within our density func- 
tional approach is generally clear, as shown by Tarazona 
et al}^. in their study of bulk systems. However, con- 
sidering such spherical structures in combination with 
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surfaces will require significant additional computational 
effort due to the further reduction of symmetry. Another 
open point is the presence of a vapor-liquid transition in 
the system. Although the present mean-field approach 
(where the anisotropic contribution to the functional can- 
cels out in a homogeneous, isotropic state) does not reveal 
any condensation transition in the bulk, our preliminary 
results from a modified mean-field approach (see Fig. [5]) 
indicate that there is, at least, a strong tendency for con- 
densation. This is in qualitative agreement to what has 
been found in another recent study^Si of Janus systems 
(where, notably, the bulk condensation transition is ac- 
companied by micellization). In the context of surface 
systems, a bulk condensation transition could have im- 
portant consequences, since it would enable a wetting 
transition (on top of the structures already observed). 
Therefore, it would be very interesting to improve the 
present density functional approach for the surface sys- 
tems beyond the mean-field level. Finally, a further in- 
teresting issue concerns the impact of curvature of the 
substrate on the self-assembled structurea^i^. Work in 
these directions is in progress. 
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contributions to the full grandcanonical density func- 
tional (|2.1ip . as well as the contribution of the exter- 
nal potential. Furthermore, /Xqx is the corresponding ex- 
cess chemical potential, A denotes the surface area, and 
Pidcai is the ideal-gas density (of a homogeneous system) . 
There are two constraints for the mixing parameter a. 
First, one must ensure in each step that (z) < 1, be- 
cause of the logarithmic term in Eq. (I2.12[) . Second, the 
convergence should be fast enough. In our calculations 
we used the following strategy (which partially follows 
that proposed in Ref . |41| ) : Within the first iteration steps 
we calculated the grand potential for several values of the 
mixing parameter. The result was fitted (using a cubic 
fit) to find the value ao where the functional becomes 
minimal. This value ao was then used to define the new 
density distribution according to Eq. (IA1[) . After few it- 
erations with this procedure (each time updating ao), we 
kept the value of ao constant in the further iterations to 
minimize numerical noise. 

As a criterion for convergence we monitored the differ- 
ence between the old and new density distribution. The 
iteration is stopped when the integral over the absolute 
values of the difference of both distributions becomes less 
than 10~^. This method yields a better accuracy in terms 
of the boundary conditions, see below. 

To calculate the weighted densities, we rewrite 
Eq. (|2.15l) in Fourier space by using the convolution the- 
orem. This yields 

n„ (z) - FFT-i (FFT (p (z)) ■ FFT {u^ (z))) , (A3) 

where FFT stands for the fast Fourier transformation^. 
The analogue of Eq. (|A3p holds for the vectorial weighted 
densities (z). In practice, the weight functions are 
multiplied with appropriate factors to allow for a 
polynomial interpolation within the integrals. 



Appendix A: Numerical implementation 

In our study, we carry out the minimization numeri- 
cally, using a one-dimensional lattice with a total length 
of = 60a and a discretization of 512-1024 points per 
sphere diameter a. We employ a simple iterative algo- 
rithm. Specifically, the density profile in step n is given 
by 



p" (z) = (1 - a) (z) -I- ap"™ (z) 



(Al) 



where a is a mixing parameter (0 < a < 1) interpolating 
between the old {p"~^) and new (p"°"') density profile. 
The new density distribution is calculated via the varia- 
tional functional derivative 



Pidoai exp 



'13 



1 SF,^[p] 



A Sp{z) 



(A2) 

In Eq. ()A2I) Fex[/o] is the excess free energy functional 
including both, the HS and the anisotropic, interaction 



There are two exact boundary conditions, which 
give us the opportunity to control the numerical ac- 
curacy. First, there is the contact theorem for planar 
walls^ 



P=-J dzp{z) J d2a;a(z,u;)^0'^"f (z,u(a;))^ 



which reduces to 



p (z = R+) = I3P 



bulk 



(A4) 



(A5) 



in the case of a hard wall at z = 0. This theorem ex- 
presses the balance of forces, that is, the force per unit 
area exerted on the wall equals the bulk pressure. The 
second condition is that the density profile far away from 
the wall approaches the bulk density (in an isotropic sys- 
tem), that is 



p{z 



Pbulk- 



(A6) 



Within our calculations, these boundaries conditions are 
fulfilled with very high accuracy for all systems at neutral 
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walls (see Sec lIIIBj) . In the presence of surface fields (see accuracy, whereas Eq. (jA6p still holds. 
Sec. IIII Ct . the contact theorem is fulfilled only with less 
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